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1.1 Introduction 

With the 1998 announcement of new evidence for muon neutrino disappear- 
ance observed by the Super-Kamiokande experiment, the more than a decade 
old atmospheric neutrino anomaly moved from a possible indication for neu- 
trino oscillations to an almost inescapable implication. In this chapter the 
evidence is reviewed, and indications are presented that the oscillations are 
probably between muon and tau neutrinos with maximal mixing. Implica- 
tions and future directions are discussed. 

The understanding of this phenomenon is now dominated by the data 
announced by the Super-Kamiokande Collaboration at Neutrino98, of which 
group the present author is a member. Much of this report dwells upon those 
results and updates to them, and so credit for this work is due to the whole 
Collaboration, listed in the Appendix, who have labored hard to bring this 
experiment to fruition and who have been ably led by Prof. Yoji Totsuka of 
the University of Tokyo. That said, this report presents personal recollections 
and opinions of the author, particularly in matters of the previous history, 
interpretation of the present situation, and future prospects for this line of 
research. 

The phenomenon of neutrino oscillations is discussed in several other 
Chapters of this volume (II and IX in particular), and to those the reader 
is directed for derivation of the expressions utilized in this Chapter and for 
understanding of the origin and implications of neutrino oscillations gen- 
erally. Model building and implications in astrophysics and cosmology are 
likewise treated elsewhere, while in the following we focus narrowly upon the 
atmospheric neutrino anomaly, its experimental explication in terms of muon 
neutrino oscillations with tau neutrinos, and the implications of those results. 



1.1.1 Atmospheric Neutrinos 

The neutrinos under discussion in this chapter arise from the decay of pions 
and other mesons, and muons, which are produced in the earth's atmosphere[0 
The atmosphere is being constantly bombarded with cosmic rays, which con- 
sist mostly of protons, but also heavy nuclei and electrons, and even neutral 
particles. The earth's magnetic field, plus other magnetic fields cut off the 
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lower energy particles from the sun and more distant sources, so that the 
mean incoming kinetic energy is around 1 GeV. Cosmic rays with lower 
energies do not cause effects which we can directly detect on earth or under- 
ground. Particles with energies in the multi-GeV range make showers in the 
roughly ten-interaction-length-thick (vertical column density) atmosphere. 
Cosmic ray collisions with air nuclei produce pions and other particles in 
abundance, which themselves further interact or decay. This competition be- 
tween interaction and decay leads to a steeper spectrum for the decay prod- 
ucts. At energies below several GeV the muons produced in charged pion 
decay, themselves decay: 

7r+ fi^ tt" ^+ e+ fi^ r>e t'p, (1-1) 

with decay lengths of L^± = 0.056 km x Et^/GcV and = 6.23 km x 
E^/GeV. Typical pion interactions lengths (roughly 150 gm/cm^) are on 
the order of a few km depending upon altitude, angle and energy, while 
muons generally come to rest before decaying or being absorbed. Moreover 
(crucially and often ignored) the energy sharing in the decays is such that the 
resulting neutrinos are also of nearly equal energy. These decay kinematics 
are of course well known, so the ratio of muon neutrinos to electron neutrinos 
can be calculated with rather good accuracy, about 5%, almost independently 
of the cosmic ray spectrum as illustrated in Figure 1.1. 
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Fig. 1.1. The calculated ratio of the fluxes of atmospheric muon neutrinos to elec- 
tron neutrinos, versus neutrino energy. Figure from HondafM]. 



Precise neutrino flux calculations (few percent) from man-made sources 
are difficult if not impossible, as indicated in other chapters in this book. 
The problem is even more difhcult for the atmospheric neutrinos, since the 
absolute magnitude of the incoming cosmic ray flux is not well known, being 
uncertain at present to perhaps 25%. |^. The calculations of the atmospheric 
neutrino flux in the few GeV energy range require not only the input cosmic 
ray flux, with appropriate modulation to account for solar cycle variation 
and geomagnetic field, but also details of nucleon-nucleon and meson-nucleon 
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interactions, not all of which have been well measured. The neutrino flux 
calculations also lead to muon flux predictions, and these can be (and have 
been) compared to data, though the appropriate data for low energy muons at 
high altitude, as recorded in balloon measurements, is sparse and imprecise. 
Typically these calculations incorporate the approximation that the incoming 
cosmic rays, the secondaries and even the neutrinos all travel in the same 
direction. This is no doubt not a serious compromise in the few GeV energy 
range, but has some effects in the energy range of a few hundred MeV. At 
this time new calculations are in progress]^, but the computer time required 
to do the full simulation is still a limiting factor and the job has yet to be 
done definitively. 

Several features of the neutrino flux are worth mentioning. The muon neu- 
trino flux can be approximated as a power law with spectral index 7 ~ —3.70 
for energies between about 10 GeV and 100 TeV. The electron neutrino (and 
anti-neutrino) fluxes which largely arise from muon decay, fall off more swiftly 
above several GeV, with strong angle dependence. As illustrated in Figure 
the ly^ to Ve flux ratio falls to a few percent at the higher energies, where 
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the i^e's are mostly produced in kaon decay, ■K°e'^Vf. (4.82% BR). See 

Figure 1.2 for atmospheric neutrino spectra of all three flavors, as expected 
under several assumptions of oscillations (two- flavor oscillations only). 

There is a significant zenith angle variation in the atmospheric neutrino 
flux, more prominent at higher energies, called the "secant theta" effect. This 
is simply due to the fact that those pions and muons which are produced by 
incoming cosmic rays with trajectories nearly tangential to the earth have 
more flight time in less dense atmosphere, and hence more chance to decay. 
Thus there is a peak, becoming is more prominent at higher energies, near the 
horizontal arrival direction in the atmospheric neutrino angular distribution. 
This peak is symmetric about the horizon for any location except at the 
lowest neutrino energies, below around 400 MeV , where geomagnetic effects 
spoil the symmetry somewhat. 

The atmospheric neutrino energies practically accessible in underground 
experiments range from the few tens of MeV to the 1 TeV range, and the 
flight distances from roughly 20 km for down-coming neutrinos to 13, 000 km 
for those traversing the earth from the far side. The neutrino cross section is 
sufficiently small that there should be negligible attenuation of these neutri- 
nos: the attenuation is roughly 2.4 x lQ~^Ey/GeV for neutrinos traversing 
the earth's core, and thus negligible for any energies below about 100 TeV . In 
consequence of the large dynamic range in both energy and flight distance, the 
atmospheric neutrinos are potentially sensitive to oscillations over a range of 
mass squared differences from about lO""* to 10~^ eV^ , as is discussed below. 



1.1.2 Initial Indications 



We will not dwell upon the past history but note that the atmospheric neu- 
trino anomaly has been around for some time, roughly since the mid-1980's. 
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Fig. 1.2. Neutrino fluxes of all three flavors (dots: v^, dashes: t/f,), long-short 
dashes Vt-) in the presence of i^t oscillations with maximal mixing and 

Arri^ = 10"^ 10"^■^ 10"^ eV^^ from left to right. Figure from Stanev||. 



Indeed the first notice of something pecuhar in the atmospheric neutrino data 
stems from the 1960's when the seminal underground experiments in South 
Africa||l^ and South lndia[0 first detected the natural neutrinos and ob- 
served somewhat of an absolute rate deficit, but not convincingly as the flux 
predictions were rough and the statistics small. 

A new round of instruments were built beginning in the late 1970's to 
search for nucleon decay as predicted by the (soon to be discarded) SU(5) 
unification model. The problem with atmospheric neutrinos, a background 
to nucleon decay searches, became serious after the activation of the first 
large water underground Cherenkov detector, the 1MB experiment, and by 
1983 the realization that the number of events containing muon decays was 
less than expected]!^. Soon this was confirmed by the second large water 
detector, the Kamioka experiment jl^ , which group extended the results with 
good particle identification giving a redundant measure of the relative muon 
deficit (as also did 1MB). Some members of the IMBlQ and the Kamioka JTst 
groups began to suggest, at least in private, that oscillations were the cause 
of the deficit, but that conclusion was not widely taken seriously for nearly 
ten years. Indeed the first published interpretations of the anomaly as due to 
oscillations were largely from outside the experimental groups To be fair 
though, it seems to be the Kamioka group who first seriously believed that 
the anomaly was due to oscillations and not simply a detector or background 
problem. 
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The deficit is characterized usuaUy as an R value, the double ratio of 
muon to electron neutrinos, observed to expected. The effect was large, the 
observed R at about 2/3 of the expected value. With the initial evidence, the 
oscillations could have been from muon neutrinos to others (eg. i^r or a new 
neutrino) or between the muon and electron neutrinos themselves. It was the 
ratio that was in deficit: one could not be sure whether there was an excess of 
electron neutrinos, a deficit of muon neutrinos, or some of both. This led to 
suggestions of other possible "physics" causes, such as nucleon decay favor- 
ing electron modes (since the anomaly was not initially detected above the 
nuc leon mass), or an excess of extraterrestrial electron neutrinos. See Table 



1.1 for a graphical summary of the situation. There were also suggestions 
of systematic problems, such as problems in muon identification, something 
wrong with flux calculations or neutrino interaction cross sections, entering 
backgrounds, or generic problems with the water Cherenkov detectors. 

Over the intervening years between the emergence of this "atmospheric 
neutrino anomaly" , as it became known, and the 1998 SuperK announcement, 
a great deal of effort went into study of these possible systematic causes of 
the anomaly. One troubling concern was that two European experiments, 
the NUSEX[0 and the Frejus|Q Detectors, did not observe any anomaly. 
Hence some people suspected a peculiarity of water as a target or with the 
employment of the Cherenkov radiation in vertex location. However, not only 
were the statistics of the European detectors relatively small, but as indicated 
by more recent work from a similar type of instrument in the US, the Soudan 
II detector the presence of a surrounding veto counter is vital for the 
more compact type of slab detectors. As well, the MACRO experiment [p2[ 
has elucidated the non-negligible production of low energy (hundred MeV) 
pions by nearby cascades in rock, which particles enter cracks in non-hermetic 
detectors and appear to be neutrino interactions. In any case the Soudan II 
with now significant exposure (4.6 kilot on-y ears) finds an R value close to 
that of SuperK (and 1MB and Kamioka)|lg. 

Figure shows these R values for the several experiments. Note that un- 
der the assumption of no oscillations all experiments should record R = 1.0, 
but if oscillations are taking place, then the R should be reduced but not 
necessarily the same value for all experiments as it depends upon the en- 
ergy range being studied. The European detectors did have higher thresholds 
which may partly explain their failure to detect the anomaly. 



1.2 The Super-Kamiokande Revolution 

We now proceed to summarize the evidence for oscillations which while de- 
pending largely upon the SuperK experiment, has important confirmation 
and consistency of results from Soudan II and MACRO (and consistency 
with previous smaller experiments as well, except as discussed below). Be- 
fore going on it may be worthwhile to point out what permitted the big 
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Fig. 1.3. The double ratio 7? of muon to electron neutrino events, data divided 
by expectations for various underground atmospheric neutrino detectors. From 
A. Mann||] 



break-through with SuperK, which may not be obvious. The increase in size 
of detector, from near kiloton fiducial volumes for Kamioka and Soudan, and 
three kilotons for 1MB to the twenty-two kilotons of SuperK is not the whole 
story. As will be seen below, the most striking progress comes from the record- 
ing of muon events with good statistics in the energy region above 1 GeV. 
This is due to detector linear dimensions as well as gross target volume: muon 
events with energy more than 1 GeV and thus 5 m range were not likely to 
be fully contained in the Kamioka detector (or the 1MB detector). SuperK 
in contrast has decent muon statistics up to almost 5 GeV, and this is turns 
out to be crucial. 

The most important data to be discussed below is the "fully contained" 
{FC) single-ring event sample, consisting of those events in which both 
the neutrino interaction vertex and resulting particle tracks remain entirely 
within the fiducial volume. For these events the relativistic charged particle 
energy and direction are well determined. We shall use the notation FG for 
the single ring events, which are about 2/3 of the total, arising mostly from 
quasi-elastic charged-current interactions in which the recoil nucleon is not 
seen. The multi-ring events have not yet been much used in analysis due to 
the ambiguous interpretation of overlapping rings from track segments in a 
Cherenkov detector (except as discussed below in the case of tests distinguish- 
ing the muon neutrino's oscillating partner). Moreover, the multi-pion final 
states are not modeled reliably in simulations as yet, and there are further 
complications of final state nuclear scattering as well. 

There are also "partially contained" (PC) events, in which a muon exits 
the fiducial volume from a contained vertex location. Such events are useful 
even though the total energy is not known, the energy observed being a 
lower limit. Of course this is the case even with FG events, though to a 
lessor degree, because the observed particles are not of the same energy (or 
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direction) as the incident neutrino, which of course is what one would like to 
know. 

The particle types are identified by pattern recognition software, now 
well tested and verified by experiment with known particle beams at the 
accelerator . Fortunately most of the contained events show single (Cherenkov 
radiating) tracks in which the identification is quite clean (at the 98% level). 



as illustrated in Figure 1.4. To be clear and cautious we usually refer to 



the reconstructed events as "muon-like" and "electron-like", though a safe 
approximation is that these represent muon and electron neutrino charged- 
current interactions. 
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Fig. 1.4. Particle identification parameter distribution of the SuperK fully con- 
tained single ring data, and Monte Carlo simulation, illustrating the electron-hke 
and muon-like separation [pi[l. 




The other two categories of events which we shall discuss are the through- 
going upwards-moving muons (UM), produced by neutrino interactions in 
the rock or outer detector, and which are coming from directions below the 
horizon (as those from above the horizon can be confused with down-going 
muons from cosmic ray interactions in the atmosphere near overhead). An- 
other category of event is the entering-stopping muon (SM). It is useful that 
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these event categories probe approximately three different energy ranges of 
neutrinos: FC ~ 1 GeV; PC and SM ~ 10 GeV; UM ~ 100 GeV. This 
is illustrated in Figure It should be understood that as far as we know, 
these neutrinos are all produced in the upper atmosphere by cosmic ray in- 
teractions, and are reasonably well described by models in content, energy, 
and angular dependence (to a few percent) 
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Fig. 1.5. Event rates as a function of neutrino energy for fully contained events 
{E < 1.3 GeV and E > 1.3 GeV), stopping muons (and similarly partially contained 
events), and through-going muons. From EngelQ. 



We shall not take up limited space here with the description of the Su- 
perK detector, which is well documented elsewhere. The interested reader 
would do well to look at some of the theses from SuperK, which are available 
on the webp5|. The short summary is that the SuperK detector consists of a 
large stainless steel cylinder (37 m high by 34 m diameter inside the inner de- 
tector) containing a structure holding 13,142 large (20 inch diameter) photo- 
multipliers . With extremely high photo-cathode coverage (40%), nearly an 
acre of photocathode and ten times more pixels than any earlier instrument, 
the instrument possesses a remarkable sensitivity of roughly eight photoelec- 
trons per MeV of deposited (Cherenkov radiating) energy. The latter permits 
detection of events down to < 5 MeV, so for the present discussion detection 
efficiency versus energy is not important because the events we are discussing 
are all above ~ 100 MeV. The inner volume is also well protected by a 2- 
m-thick, fully-enclosing veto Cherenkov counter, populated by 1800 recycled 
1MB (8 inch) photo- multipliers with wavelength shifting collars. The inner 
"fiducial" volume is further taken as 2 m inside the inner photo-multiplier 
surface, resulting in the 22.5 kiloton volume used for most reported data. 
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The SuperK oscillations claim was first formally presented to the physics 
community in June 1998 at the NEUTRIN098 meeting in Takayama, Japan. 
The data were presented in several papers to the community j2^,^,^, build- 
ing upon past data from Kamioka[|l5j and IMB|l4j, and culminating in the 
claim of observation of oscillations of muon neutrinos, published in Physical 
Review Letters in August 1998 We now proceed to review the evidence, 
which has changed little except for new indications that the i^^ oscillating 
partner is probably the Vt, and not a hypothetical sterile neutrino. 



1.2.1 Up-Down Asymmetry 

One way to look at the FC (and PC) data is in terms of a dimensionless up- 
to-down ratio, difference over sum (which has symmetrical errors in contrast 
to just up/down) |]3l[. Downwards going neutrinos have flown ^ 20 — 700 km 
while up going neutrinos have traveled ^ 700 — 10,000 km. The angle be- 
tween the neutrino and observed charged lepton is on average of the order of 
40°/ ^jE^^jGeV and the typical observed energy is a half the neutrino energy. 
Thus the mixing of the hemispheres of origin of the events is important only 
for the lowest energies (below roughly 400 AleV). This asymmetry quan tity 



is exhibited as a function of charged particle momentum in Figure L6, for 
both electrons and muons, with the PC data shown as well (for which we 
know only a minimum momentum), from an exposure of 70.4 kiloton years in 
SuperK. One sees that the electron data fit satisfactorily to no asymmetry, 
whilst the muon data show strong momentum dependence, starting from no 
asymmetry and dropping to about -1/3 above 1.3 CeV. 

From this Figure alone, without need for complex and often opaque Monte 
Carlo simulations, assuming the cause of deviation from uniformity to be due 
to neutrino oscillations, one can deduce that: 

1. The source of the atmospheric neutrino anomaly is largely due to disap- 
pearing muons, not excess electrons. 

2. There is little or no coupling of the muon neutrino to the electron neutrino 
in this energy/distance range. 

3. The oscillations of the muon neutrinos must be nearly maximal for the 
asymmetry to approach one third. 

4. The scale of oscillations must be of the order of 1 GeV^/200 km, plus or 
minus a factor of several. 

In fact, as seen by the dashed lines overlying the data points, the simu- 
lations do produce an excellent fit to the muon neutrino oscillation hypoth- 
esis, while the no-oscillations hypothesis is strongly rejected. The latter is 
so strong that statistical fluctuations as the cause of the deviation are com- 
pletely improbable; one must look for systematic problems in order to escape 
the oscillations explanation. 

One concern for some people has been the fact that the asymmetry is 
indeed maximal, which makes it appear that we are very lucky that the earth 
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Fig. 1.6. The up-to-down asymmetry for muon (2486) and electron (2531) single 
ring fully contained and partially contained (665) events in SuperK, from 1144.4 
days of live time (analyzed by 6/00), as a function of observed charged particle mo- 
mentum. The muon data include a point for the partially contained events (PC) 
with more than about 1 GeV. The hatched region indicates no-oscillation expec- 
tations, and the dashed line — Ut oscillations with Z\m^ = 3.2 x 10~^eV'^ and 
maximal mixing |2ll]. 



size and cosmic ray energies are "just so" to produce this dramatic effect. 
This appears to this author to fall in the category of lucky coincidences, 
such as the angular diameter of the moon and sun being the same as seen 
from earth. There is another oscillations related peculiar coincidence that 
the matter oscillation scale turns out to be close to one earth diameter, and 
this depends upon the Fermi constant and the electron column density of the 
earth. The phase space for "coincidences" is very large, and we humans are 
great recognizers of such patterns. 

1.2.2 Neutrino Flux Dependence Upon Terrestrial Magnetic 
Field 

The effect of the earth's magnetic field on the atmospheric neutrino flux is 
a little complicated, but only important for very low energies. For example 
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for energies of a few GeV, the magnetic field provides some shielding from 
straight downwards going charged cosmic rays in regions near the magnetic 
equator. For higher energies and incoming trajectories near the horizon, the 
magnetic field still prevents some arrival paths. As the SuperK detector loca- 
tion is not on the magnetic equator the effect is not up-down symmetric, and 
this spoils the symmetry otherwise expected from the neutrinos about the 
horizontal plane. However, the effects are mostly limited to neutrino ener- 
gies below about 1 GeV, corresponding to cosmic ray primaries below about 
10 GeV. The picture is made a bit more complicated by the earth's mag- 
netic field not being a nice symmetrical dipole. Fortunately there are good 
models of the magnetic field, and the people who have made flux calcula- 
tions take these effects into account, though (in the past) largely through a 
simple cutoff momentum versus location. More recent calculations trace par- 
ticles backwards in the magnetic field and determine trajectories that escape 
to infinity |p3| . Lipari has, however, recently shown that the double humped 
cosmic ray spectra seen in the AMS experiment, with a space borne mag- 
netic spectrometer in low earth orbit, may be due to particles in trapped 
orbits[^. Moreover, Lipari points out that there are hints in the AMS data 
of a North-South asymmetry, which could bias the neutrino flux calculations, 
and even pull the derived value of Am? . However, it should be emphasized 
that the effect of such variation from simple expectations will only bias the 
lowest energy data in SuperK (roughly below 400 MeV), and analysis has 
demonstrated that the results quoted herein are stable against raising the 
acceptance energy for the data sample. 

The SuperK group has published a paper examining the azimuthal 
variation of the SuperK data (±30 deg about the horizon) for intermediate to 
higher energies (400 — 3000 MeV) , in an energy region where the calculations 
are thought to be reliable. Indeed the SuperK data do exhibit significant 
variation from uniformity while fitting the flux predictions very well, giving 
one confidence in the modeling [2^. 

1.2.3 Natural Parameters for Oscillations: L/E 

In an ideal world, one would assuredly study these data as a function of dis- 
tance divided by energy, L/E, since that is the parameter in which one ex- 
pects to see oscillatory behavior. For two-neutrino mixing with mass squared 
difference, Am? , and mixing angle 9, the probability of a muon neutrino of 
energy E^, remaining a muon neutrino at distance L is given by []63| 



However, since we observe only the secondary charged particle's energy and 
direction, badly smeared at the energies available {L/E^ smeared by about 
a factor of two), plots in which one would wish for visible oscillations can 
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at best show a smooth shdc from the no-oscillations region to the oscillating 
regime. This is illustrated in Figures |l.7| and 1.8, where the numbers of events, 



and the ratios of those numbers of events observed to those expected with 
no-oscillations, are plotted versus "L/i?"|3^, for muon and electron (type) 
events. The updated data are preliminary from the SuperK 1144 day sample. 
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Fig. 1.7. The numbers of SuperK events observed compared to predicted as a 
function of the natural oscillations parameter, L/E, distance divided by energy. The 
results are not normalized. The two peaks correspond to generally down-coming 
(left) and up-coming (right). One sees that the muon deficit begins even in the 
upper hemisphere. The shaded area indicates no-oscillations expectations, and the 
dashed line in the fit for i/^ — Vr oscillations with maximal mixing and Z\m^ = 
0.0032 eV'^M] 



The plot is not "normalized" , and we see somewhat of an excess of electron 
type events overall (-1-8%). This is a little worrisome, but acceptable since (as 
already noted) the absolute flux is uncertain to a larger value. In contrast to 
the electron data, the muon points fall relative to no-oscillations expectations 
with increasing L/E beyond about 50 km/GeV , reaching a plateau at about 
one half their initial value, consistent with maximal mixing. Muon (to tau) 
neutrino oscillations in the Monte Carlo simulation are indicated by a dotted 
lines, and fit the data reasonably well. 
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As noted, these data do not (and could not) show oscillations, due to 
convolutions washing out the oscillatory behavior. It was this smooth fall, 
however, that caused the author and some colleagues to wonder if another 
model might fit the data, one in which one component of the muon neutrino 
decays rather than oscillates with distance. Two papers suggested neu- 
trino decay to explain the atmospheric neutrino anomaly. I will not discuss 
details here, but note that in order to construct a viable model we had to 
push on all available limits, and invoke neutrino mass and mixing in any case. 
Consequently such models do not pass the economy test of Occam's Razor, 
though most annoyingly they remain not ruled out as yet. 
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Fig. 1.8. The ratio of numbers of events observed compared to predicted as a 
function of the natural oscillations parameter, distance divided by energy. The 
results are not normalized and overall there is a slight excess (about 8% compared 
to a systematic uncertainty of 25%) compared to expectations. Electrons show no 
evidence for oscillations, while muons exhibit a strong drop with L/E. This is 
consistent with — Ut oscillations with maximal mixing and Am^ — 0.0032 eV^, 
as indicated by the dashed lines from the simulation. pl[| 



One may note that detecting multiple oscillation peaks is not ruled out 
in principle for detectors such as SuperK or Soudan II, employing the at- 
mospheric neutrinos. It is a matter of recording the final state of the muon 
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neutrino charged eurrent events, including nuclear recoil, with sufficient ac- 
curacy and statistics. Detectors such as a liquid argon device of the ICARUS 
type are claimed to have the resolution, if large enough. Soudan II has ap- 
parently good enough resolution to accumulate a "golden sample" in which 
the nuclear recoil is detected, permitting reconstruction of incident neutrino 
energy and direction. Unfortunately Soudan II does not have enough mass 
to achieve definitive statistics in a practical observing period||l^. Another 
possibility is that SuperK with enough exposure and more highly developed 
analysis would be able to accumulate an adequate sample of events in which 
the recoil proton is detected above Cherenkov threshold. At the moment none 
of the above promise success. 

Considering future experiments, the attempt to discern oscillations as a 
function of L/ E is one area in which improvement may indeed be made. 
The MINOS|62l detector in Minnesota with a neutrino beam from Fermilab, 
and the large detectors to be constructed in Gran Sasso, ICANOE|^ and 
OPERA []62[, detecting a neutrino beam from CERN, give some hope of being 
able to yield oscillatory plots. A hypothetical detector, such as a megaton 
version of the Aqua-RICH instrument studied by Ypsilantis and colleagues 
could have the resolution to see a multi-peaked L/E plot|3|,^. Nearer to 
technical development, the proposed MONOLITH experiment would consist 
of a 30 kiloton pile of magnetized iron and tracking detector layers, and 
employ the cosmic rays to detect the first dip in L /E in the upcoming muon 
flux through the earth p7[. 



1.2.4 Energy and Angle Variation 

The SuperK Collaboration's preferred method of fitting the ensemble (single 
ring) EC and PC data is to employ a test to numbers of events binned 
by particle type, angle, and energy, a total of 70 bins. The bin choices may 
seem a bit peculiar, but they have historical precedent (they are as employed 
for Kamiokande) and though not optimal for the new data set, this choice 
permits avoidance of any statistical (or confidence) penalty for choosing ar- 
bitrary bins. The fit employs a set of parameters to account for potential 
systematic biases. Details cannot be presented here, but it has been shown 
that the numerical results are quite insensitive to the selection of the param- 
eters or their supposed "errors" (except for the overall normalization) ||38| . 
This method of systematic error handling has been shown to be equivalent 
to employment of the correlation matrix of parameters |32|. 



Figure illustrates the data plotted for two energy intervals {sub — GeV 
and multi — GeV, less or more than 1.3 GeV) for single track events identified 
as either electron-like or muon-like. The partially contained data is displayed 
with the multi-GeV muon data. The data are shown as a function of the 
cosine of the zenith angle, with -1-1 being down-going. One sees that the data 
very well fit the curves from the Monte Carlo simulation, at the values gotten 
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from the grand ensemble fit, Am^ = 0.003 eV^ and si-n?{26) = 1.0. The limit 
on Am^ is 0.002 to 0.007 eV'^ , and sin2(26') > 0.85 at 90% confidence level. 
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Fig. 1.9. Cosine of zenith angle distributions of the contained and partially con- 
tained event data for two different energy ranges (above and below 1.3 GeV), 
electron and muon single ring events. 1144 live days of SuperK data (preliminary 
analysis) are indicated by dots with statistical error bars. The solid line shows 
no-oscillations simulation result, and the hatched line that for oscillations between 
muon and tau neutrinos with best fit Am? = 0.0032 eV^ and maximal mixing pl[]. 



The results of the fits are often presented in terms of an inclusion plot, 
showing an acceptable region(s) in the space of mixing angle (sin^ 29) and 
mass squared difference [Am?), as presented in Figure 1.10. The Am? value at 

minimum has moved a little upwards with accumulated statistics, though 
not much, (good news for long baseline experiments anyway) but remains 
uncertain to about a factor of two. 

It is noteworthy that the earlier indications of and constraints upon oscil- 
lations parameters from Kamiokande, 1MB and Soudan gave somewhat larger 
values of Am? . All of these results depended upon fitting the R value, since 
no angular distribution was discerned (due to limited statistics and lower 
mean energy due to containment). Later Kamiokande data did show angular 
variation in the PC data, but not statistically compellingly. For some reason 



16 John G. Learned 




68% C.L 
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Fig. 1.10. Inclusion plot, showing the regions for 3 levels of statistical acceptability 
in the plane of mixing angle and mass squared difference for Vfi — u-r oscillations. This 
is from SuperK contained and partially contained event 1144 day data preliminary 
analysis of June 2000 ||l[. 

not fully understood, the fits using R alone all seem to yield higher values of 
Am? . If one has some deficit in muons without angular determination, then 
one can fit that suppression with any Am? above some threshold value by 
choosing an appropriate mixing angle. Thus the R constraints are open ended 
upwards in Am^ . Perhaps there is a systematic problem here due to predicted 
neutrino spectra, or perhaps there is some physics yet to be elucidated. This 
is not to suggest that it seems possible for the preferred two-neutrino solution 
to move much, but that more complex small effects at the < 10% level could 
be superposed on the present simple solution. Accelerator based experiments 
should clarify this issue. 

1.2.5 Muon Decay Events 

It is not often emphasized, but the original indication of the anomaly, a 
deficit in stopped-muon decays (~ 2.2 /isec after the initial neutrino event), 
remains with us, and constitutes a nice alternate sample, almost independent 
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and with quite different systematics. It is not so clean a sample (there are 
muon decays from pions produced in electron CC and all flavor NC events) 
and the statistics are lower, but the complete consistency of the muon decay 
fraction remains a reassuring complement to the energy and angle analysis 
employing track identification. 



1.2.6 Through-Going and Entering-Stopping Muons 

Another cross check comes from the UM and SM samples, which are par- 
ticularly attractive because the source energies are factors of 10 and 100 
higher and the detector systematics rather different (for example, the tar- 
get is mostly rock not water). A drawback to these samples is that one is 
restricted to using muons arriving from below the horizon, due to the over- 
whelming number of down-going cosmic ray muons penetrating the mountain 
(at 50,000 times the rate in SuperK). 

In going from the earlier instruments to SuperK, however, the gain is not 
so great (the 1200 of SuperK being about a factor of three more than 
the previously largest underground instrument, 1MB, for example), since the 
rate of collection of through-going muons depends upon area not volume. 
However, the much greater thickness of the detector (and the efficient tagging 
of entering and exiting events in the veto layer) yields many more entering- 
stopping (SM) events. 

The flux angular distribution derived from 1260 UM events from below 
the horiz on, ea ch with more than 7 m track length in the detector, is shown 



in Figure |1 . 1 1| , where one sees that the angular distribution is nicely consis- 
tent with oscillations and not with no-oscillations. However, since much of 
the effect is close to the horizon, where oscillations for the energies in ques- 
tion are just setting in, one worries about contamination of the near horizon 
events with in-scattered events from the much greater numbers of down going 
muons. There is no room for detail here, but SuperK does perform a small 
background subtraction (9 events of 247, or 3.6% in that one bin) for events 
within 3 degrees of the horizon, but otherwise finds no evidence for significant 
contamination |]30| . 

In SuperK the SM sample was predicted to be 33-42% of the UM sample. 



as indicated in Figure 1.12 , yet in fact SuperK sees only about 24% ± 2%. 
Fitting the data to the oscillation hypothesis, one can make the now usual 
inclusion plot, which shows that the UM and SM results are completely in 



accord with those from the FC and PC data, see Figure 1.23. However, as the 
statistics are smaller and the physics leverage not as great, the muon result 
does not add much to the FC and PC constraints, though it does stiffen 
the lower bound in Am? . The joint fit to the UM and SM data alone yields 
a /ndf — 35.4/15 and 13/13 for the cases of no-oscillations and — Vr 
oscillations. 

One may note that earlier experiments, such as 1MB, with a final sample 
of 647 events, no veto counter and less mature flux calculations, did not 
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Fig. 1.11. The prehminary flux calculated from 1260 upcoming muons during 
1138 live days of SuperK data, with 9 events background subtracted in bin near- 
est the horizon, cos(6') — 0. Error bars show statistical errors only. Expecta- 
tions for no-oscillations (solid line) and the best fit for — Vr oscillations with 
Arn^ = 0.0032 eV^ and maximal mixing (dashed line), are shown pl|]. 

find any net deficit in the UM sample, nor significant deviation from no- 
oscillations expectations. A similar case obtained with other smaller data 
sets. Indeed, one may note that on the strength of the SuperK UM data 
angular distribution alone, one would hardly be making discovery claims. All 
UM data from 1MB and Kamioka are and were in accord with the present 
results, but not demanding of the oscillations conclusion. 

There is a lengthy tale about an SM/UM analysis from the 1MB experiment p^, 
which claimed an exclusion region very close to the now preferred solution. 
The 1MB stopping muon sample was small and it was not clean due to lack of 
a veto layer. More importantly the interpretation seems to have been flawed 
due to older flux models and Monte Carlo simulations. Work is in progress 
to reassess the old data with new flux calculations and an updated quark 
model p9|. Thus there remains a cloud upon the horizon, but one which may 
fade away upon reanalysis. It might also be worth recalling that in the 1980's 
the absolute rate of upcoming muons as measured in the 1MB, Kamioka and 
other detectors, agreed with calculated rates employing then available flux 
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Fig. 1.12. The flux of stopping muons versus zenith angle in 1117 days of SuperK 
data. There are 311 events and a background of 21.4, with track length greater than 
7 m, or about 1.6 GeV. The large overall deficit in the data (dots with statistical 
error bars) compared to the no-oscillations (solid line), is less significant than it 
appears because of the 20% uncertainty in absolute fluxEl|]. 



calculations. Also at that time, peculiar angular distributions which fit no 
expectation were reported at conferences from the MACRO and Baksan de- 
tectors. These results all tended to give pause to claiming oscillations as the 
resolution of the atmospheric neutrino anomaly. These concerns were swept 
away by the clean and statistically convincing FC muon angular distributions 
from SuperK. 



1.2.7 The Muon Neutrino's Oscillation Partner 

Given that the muon neutrino is oscillating, is it oscillating with a tau neu- 
trino or a new sterile neutrino which does not participate in either the charged 
(CC) or neutral current (NC) weak interaction? Fortunately there exist sev- 
eral means to explore this with SuperK data. The NC interactions should 
show an up-down asymmetry for sterile neutrinos but not for tau neutrinos 
(since the NC interactions for all ordinary neutrinos are the same). Another 
avenue for discrimination is that sterile neutrinos would have an additional 
oscillation effect due to "matter effects" . The consequence would be a unique 
signature in the angular distribution of intermediate energy muons, as illus- 



trated in Figure 1.13 



Early SuperK efforts focused upon the attempt to collect a clean sample of 
7r° events. As it turns out, this was frustrated because the rings (from the two 
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S> 0.8 




Fig. 1.13. Survival probability P(i'^i I'fi) as a function of the zenith angle in 
the cases of maximal mixing of V/j, with (upper panel) , (middle panel) and i/s 
(lower panel). For |Zim^| — 5- 10"'' eV^ the curves correspond to neutrino energies 
20, 40, 60 and 80 GeV. The dashed curves are calculated with the approximation 
of constant average densities in the mantle and in the core of the earth. From 
Lipariliol. 
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decay 7s) cannot be separated at energies above ~ 1 GeV, and in net there 
are not so many reconstructed events as to permit a good discrimination. In 
fact the absolute rate is consistent with expectations, but the cross section 
is uncertain to about 20% making the hint at tau coupling not significant. 
The K2K experiment should soon measure this cross section to perhaps 5% 
however, making the n° rate a useful discriminant. 

More recently, tests have been devised employing a multi-ring sample 
(MR), the PC event sample, and the UM sample, all of which are inde- 
pendent of the single ring FC sample which yields the strongest oscillation 
parameter bounds. 

The MR sample is cut by energy (> 1.5 GeV) and the requirement of the 
dominant ring being electron-like to enhance the NC content of the sample. 
A test parameter is constructed from the ratio of events fr om w ithin sixty 
degrees of the zenith and nadir. This is illustrated in Figure 1.14 , where one 
sees consistency of Vt and disfavoring of v sterile ■ 




Fig. 1.14. (a) Zenith angle distributions of the multi-ring events satisfying cuts as 
described in the text. coaO is -1-1 for down-going events. The black dots indicate the 
data and statistical errors. The solid line indicates the prediction for tau neutrinos, 
and the dashed for sterile neutrinos with (Am^ , sin^26) = (3.3 x W~^eV^ ,1.0). 
These two predictions are normalized by a common factor so that the number of 
the observed events and the predicted number of events for i/f^ ^ Vr are identical, 
(b) Expected up/down ratio as a function of Arr? . Horizontal lines indicates data 
(solid) with statistical errors (dashed). Black dots indicates the prediction for tau 
neutrinos, and the empty squares for sterile neutrinos, both for the case of maximal 
mixing [klf. 
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The PC sample can be cut on energy (requiring > 4 GeV) in order 
to achieve a higher neutrino source energy, and the upwards going number 
compared to downwards number of events. In this instance one is seeking 
matter effects, and the results are shown in Figure 1.15| , indicating again a 
preference for Vt over VsterUe- 




COS© Am (eV ) 

Fig. 1.15. (a) Zenith angle distributions of partially contained events satisfying cuts 
of Evia > 5 GeV. cos© is +1 for down-going events. The black dots indicate the 
data and statistical errors. The solid line indicates the prediction for tau neutrinos, 
and the dashed for sterile neutrinos with (Arn^ , sin^26) = (3.3 x 10~"^eV'^, 1.0) 
(b) Expected up/down ratio as a function of Am^. Horizontal lines indicate data 
(solid) with statistical errors (dashed). Black dots indicates the prediction for tau 
neutrinos, and the empty squares for sterile neutrinos, both for the case of maximal 
mixing Pi . 



For the muons, a near horizontal number can be compared to a number of 
nearly straight upcoming events for another test of matter oscillations. This 
is presented in Figure 1.16. 

Finally the three tests are combined in a single t^st for the case of 



Vt and the two cases for 



^sterile 

heavier or lighter than i-^. 



This is 



presented in Figure 1.17, where one sees that the entire region in mixing 
parameter space is eliminated for sterile neutrinos at more than the 99% 
confidence level, whilst the Vj- case fits perfectly [p] 

As to muon neutrinos coupling to electron neutrinos, SuperK can say 
only that the Am^ is out of range on the low side or that the sine of the 
mixing angle is less than about 0.1, if Am^ is large. As indicated in the 
earlier plots with up-down asymmetry, there is surely not much mixing in this 
energy range. One clever scenario 1 45 has the electron neutrino oscillation loss 
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COS© Am (eV ) 

Fig. 1.16. (a) Zenith angle distributions of the upward moving through-going 
muons. COS& is -1 for vertical up-going events. The black dots indicate the data 
and statistical errors. The solid line indicates the prediction for tan neutrinos, and 
the dashed for sterile neutrinos with (Z\m^sm^26') = (3.3 x 10"^eV"^ 1) (b) Ex- 
pected up/down ratio as a function of Am^ . Horizontal lines indicates data (solid) 
with statistical errors (dashed). Black dots indicates the prediction for tau neu- 
trinos, and the empty squares for sterile neutrinos, both for the case of maximal 
mixing 

being just compensated by muon neutrinos splitting their oscillations between 
electrons and taus. This seems to be ruled out by higher energy SuperK data 
however [p2| . Further discussions of (many) other scenarios of oscillations can 
be foimd in the Chapter IX. 

1.2.8 Sub-Dominant Oscillations 

For all of the foregoing we have considered only two flavor oscillations. A 
simple three flavor analysis can be made employing the assumption that solar 
neutrino oscillations are driven by a much smaller mass difference than found 
for — Vt, mi < TO2 << m^. Then, using the standard notation for three 
neutrino MNS mixing matrix (see Chapter IX), the oscillation probabilities 
can be written as: 



p 


= l-4sin2(26i23)cos2(6ii3)(l-W(6l23)cos2( 


613)) sin^ {1. 27 Am'^L/E),{l -5) 


P 


= cos*(6ii3) sin2(26l23) iA^[\.TI Arn^LlE), 


(1.4) 


p 


= sin2(26ii3) sin2(6'23) %\^[\.T1 Am^LlE), 


(1.5) 


p 


= sin2(26'i3) cos2(6l23) ^^^{l.TI Ani-LjE), 


(1.6) 
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Fig. 1.17. Excluded regions for tiiree alternative oscillation modes, (a) ^ Vt, 
the light (dark) shaded region is excluded at 90(99)% C.L., (b) i/^ ^ i^sterUe with 
Am^ > 0, the whole region shown in this figure is excluded at the 99% C.L., and 
(c) Arn^ < 0, whole region is excluded at 90% C.L., the dark shaded region is 
excluded at more than the 99% C.L.. The thin dotted (solid) line indicates the 90 
(99)% C.L. allowed regions from the FC single-ring events [}fl[. 
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P,.,^, = 1 - sin2(26'i3) nin^ [1.27 Arn^L/E), 



(1.7) 



where Am? — m^^ — m2^, and we neglect any consideration of CP or 
CPT violation. 

The SuperK FC and PC data 990 day data has been fitted with these 
equations p6[ , and the limits on sin^(26'i3) are less than 0.25 at 90% C.L., with 
Arn^ = 0.003 eV'^. The best fit value lies at sin'^{ei^) = 0.03, 5^2(6*23) = 
0.63. At this Am^ the CHOOZ resuhs allow either 5^2(6*13) < 0.03 or > 0.97 
so the latter is eliminated by the SuperK results, as illustrated in Figure 1.18 
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Fig. 1.18. Three neutrino fits to the SuperK data showing 90% and 99% al- 
lowed regions. The shaded region is the 90% exclusion region from the CHOOZ 
experiment Bf| . See text for qualifications. 



1.2.9 Non-Standard Oscillations 

As discussed in Chapter IX and elsewhere there are models of neutrino 
oscillations which result from gravitational splitting, violations of Lorentz 
invariance, and so on, which result in oscillation with a phase proportional 
to, say, L X E instead of L/E, or even with no energy dependence. The 
SuperK group has presented a fit to the data as a function of 

= sm\20) sm\pLE" ) (1.8) 



where n, /3 and sin^(20) (0.7 to 1.3) are varied. As indicated in Figure 1.19 , 
the exotic solutions are strongly disfavored, while the normal function is 
perfectly acceptable, with n = —1.06 ± 0.14. 
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Fig. 1.19. The variation of with index n in oscillating phase for muon neutrinos, 
as in const * L * E" . n = -1 corresponds to ordinary oscillations. One sees that non- 
standard solutions with n of and 1 are strongly disfavored 



1.2.10 Hypotheses to Explain Anomaly 



We conclude with a summary, presented in Table IT, of all hypotheses put 
forth to explain the atmospheric neutrino anomaly. Space does not permit 
a full discussion here, but it is the case that the SuperK data now have 
eliminated almost all alternate hypotheses to explain the results. 

The atmospheric neutrino flux calculations cannot be producing the anomaly 
since we see the effect otherwise unexplained in the muon zenith angle dis- 
tribution. The cross sections cannot produce such a geographical effect, even 
if one could find some lepton universality breaking phenomenon. Particle 
identification has been heavily studied and verified at the accelerator with 
a 1000-ton Cherenkov detector tankQ. Entering backgrounds, which would 
produce effects clustering near the outer walls, show no evidence for contribu- 
tion to the anomaly in the SuperK data. Detector asymmetrical response is 
ruled out by the observed symmetry of the electron data, as well as detailed 
calibration studies with isotropic sources. 

Extra-terrestrial neutrinos cannot be the culprits (unfortunately), since 
we see that it is a deficit of muon neutrinos causing the anomaly, not an excess 
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of electron neutrinos. Proton decay is ruled out similarly due to geographical 
dependence of the deficit, and the extent of the anomaly to too high an en- 
ergy (ruling out neutron-anti- neutron oscillations as the source as well) . The 
decay of a muon neutrino component has been discussed, and while seem- 
ingly unlikely is not totally ruled out as yet. Anomalous absorption of muon 
neutrinos by the earth, correlated with an exponential in the column density 
through the earth, is ruled out as well. This follows because the anomaly is 
not dominated by that small part of the solid angle going through the earth's 
core, whereas we see the muon deficit starting even above the horizon. Elec- 
tron neutrino mixing is ruled out as the dominant effect, again due to the 
muon angular distribution. Sterile neutrinos do not fit the data as discussed 
above. As discussed in the previous section, non-standard oscillations are also 
ruled out. 

The only hypothesis which survives, and which fits all the evidence, is 
that muon neutrinos maximally mix with tau neutrinos with a Am? in the 
range of 2 — 7 x 10^"^ eV^. It is noteworthy to this author that in all the tests 
made on the data sample, there appears to be great stability in the results 
against variations of all the parameters explored. 



1.2.11 Results from Soudan II 

The Soudan 2 detector located in an old iron mine in Minnesota, USA, con- 
sists of a vertical slab, fine-grained tracking calorimeter of 963 tons total mass. 
The cavern has a surrounding layer of proportional tubes, 2 or 3 layers thick 
on all sides. Data have been reported from 4.6 kiloton years exposure|20|. 



as illustrated in Figure 1.20. The contained events plotted are selected for 
lepton energy > 700 M eV with no visible nuclear recoil, or visible energy 

> 700 MeV and summed momentum > 450 MeV/c and lepton momentum 

> 250 MeV/c. The energy resolution is of order 20%, and the angular res- 
olution of order 20-30 degrees. In the figure the predicted number of events 
has been normalized to the electron total. One can see that with only of the 
order of 100 events of each type, the statistical significance is not great, but 
the depletion of upcoming muon events is evident. The fits to oscillations 



parameters are included below, in Figure 1.23 



1.2.12 Results from MACRO 

The MACRO detector, built primarily to seek monopoles, possesses a signif- 
icant capability, with effective mass of 5.3 kilotons, to detect through-going 
and stopping muons as well as contained and partially contained neutrino 
interactions ]4^ . The instrument, located in the deep underground Gran Sasso 
National Laboratory in Italy, consists of horizontal planes of a tracking in- 
strument. 



Figure 1.21 shows the results of analysis of partially contained data, for 



which up and down cannot be distinguished, but which shows a clear deficit 
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Table 1.1. List of hypotheses invoked to possibly explain the atmospheric neutrino 
anomaly. The first 3 columns are criteria available prior to SuperK, and the last 
4 after the 1998 SuperK publication. The hypotheses divide into 5 systematics 
issues and 8 potential physics explanations. As indicated in the text, the only 
remaining likely hypothesis is the oscillation between muon and tau neutrinos. The 
"x" schematically indicates which evidence rules out the hypothesis in that row. 
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Fig. 1.20. (left): Distributions in cos 6^ for the and flavor event samples. Data 
(crosses) are compared to the null oscillation MC (dashed histogram) where the MC 
has been rate normalized to data, {middle): Distributions of \og{L / E^) for i^e 
and i/fi charged current events compared to the neutrino MC with no oscillations; 
the MC has been rate-normalized to the Ue data, (right): Comparison of L/E^ 
distribution for data (crosses) and expectations from neutrino oscillations for 
four Arri^ values, with sin2(26') = 1.0. From Mann|20l 



compared to expectations with no-oscillations. The acceptance of such a pla- 
nar instrument is small near the horizon so most of the effect is from the 
nearly vertical events. The In-Up sample has 116 events and exhibits a de- 
pletion of 0.57 ± 0.16 compared to expectations. 
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Fig. 1.21. Distributions in cosS for MACRO partially contained events. The data 
(solid circles) are seen to fall below the null oscillation expectation in every bin of 
both samples. The dashed line shows expectations for maximal mixing and Arn^ — 
0.0025 eV^ From SurdoH. 



30 John G. Learned 



The MACRO detector also has a significant sample of upwards through- 
going muons, as illustrated in Figure 1.22| . Note that the muon energy thresh- 
old is low, in the hundred MeV range (varying with angle and entry location) . 
The overall depletion (data/expectations) is 0.74 ± 0.03 ± 0.04 ± 0.12, where 
the last term reflects uncertainty in flux and cross section[Q. While the flt to 
oscillations appears not to be perfect and the minimum in ^ lies in the un- 
phys ical region, the confldence level boundaries shown on the right in Figure 
|l.22 indicate consistency with the SuperK results. 




Fig. 1.22. a) (left): The angular distribution of upward through-going muons 
observed in MACRO. The data distribution (solid circles) differs from the null os- 
cillation expectation (shaded band) in shape and in rate, b) (right): The neutrino 
oscillation allowed regions obtained by MACRO from the upward through-going 
muons. Confidence level and experimental sensitivity boundaries are calculated us- 
ing the Feldman- Cousins method. From ^jf . 



1.2.13 Combined Evidence 



In Figure 1.23 on the left is the combined flt of FC and PC (848 live days) 
plus UM (923 live days) plus SM (902 live days) from SuperK, with results 
as indicated, constraining the oscillation parameters to be roughly 0.002 < 
AnP'leY'^ < 0.006 and 0.85 < sm^{2e) < 1.0, with the fit minimum in 
the physical region at maximal mixing and 0.0035 eV^. The boundaries for 
MACRO and Soudan 2 are shown overlying the SuperK results in the right 
hand panel. One can sec that the results of the largest instruments now 
reporting atmospheric neutrino data, water Cherenkov and two dissimilar 
tracking calorimeters, all agree on muon neutrino disappearance, and are 
consistent with oscillation between muon and tau neutrinos. I cannot do 
better than to quote Tony Mann from his 1999 Lepton-Photon conference 
plenary talk|2^: "I propose to you that congratulations are in order for the 
researchers of Kamiokande and of Super-K and more generally, for the non- 



accelerator underground physics community. For Fig. 1.23 b, Ladies and Gen- 
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tlemen, is the portrait of a Discovery - the discovery of neutrino oscillations 
with two-state mixing." 
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Fig. 1.23. a) (left): Allowed regions obtained by Super-K based upon fitting to 
FC and PC single ring events, plus upward stopping muons, plus upward through- 
going muons. b) (right): Oscillation parameter allowed regions from Kamiokande 
(thin-line boundary), Super-K (thick-line boundary), MACRO (dashed bound- 
aries), and Soudan 2 (dotted and dot-dashed boundaries). From MannpO|. 



1.2.14 Long Baseline Results 

The K2K experiment has been in operation for over one year at this writing. 
The 12 GeV KEK proton synchrotron has been equipped with a neutrino 
line, and a double detector has been built at about 100 M range, on the 
KEK campus. Neutrinos are aimed at the SuperK detector at a distance of 
250 km, and events have been recorded ||5^. Note that a 1 GeV neutrino 
would be near first minimum after 250 km if the Am^ = 0.005 eV^. 

At the time of this writing preliminary results are available, and it has 
been reported in conferences that the rate is low, consistent with the SuperK 
results, and inconsistent with no-oscillations at the 2 a level [|59[|. From the 
period of June 1999 through March 2000, the group received about 1.7 x 10^^ 
protons on target (17% of proposal request), during the equivalent of about 
five months running. The events were easily extracted from a GPS synchro- 
nized 1.5 usee time window relative to the 12 GeV KEK proton synchrotron 
fast beam spill (2.2 sec repetition period and about 5 x 10^^ protons on target 
per spill). During this run 17 associated events were collected in the SuperK 
22.5 kiloton fiducial volume, when 29.2±3.1 were expected for no oscillations. 
The expected numbers are 19.3, 12.9 and 10.9 for 3, 5, and 7 xlO^^ eV'^. Of 
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these SuperK events, 10 are single ring, one of which is identified as electron- 
like (about equal to expectations). 

Given the low energy of the neutrino beam (~ 1 GeV), and frustratingly 
low data rate at SuperK probably not much more can be expected from this 
experiment than a simple confirmation of the SuperK atmospheric neutrino 
results. The full proposal-run would collect 174 or so events with no oscilla- 
tions, and assuming oscillations, about half that. Of those, about half should 
be single ring muons, leaving perhaps 40-50 events from which to deduce the 
arriving neutrino spectrum. Given that the accelerator delivers beam about 
half a year each calendar year, one can see that it will take several years to 
achieve definitive results. 

1.3 Implications 

The ramifications of the explication of the atmospheric neutrino anomaly 
in terms of neutrino mixing and thus neutrino mass, are great and span 
the known realms of fundamental physics from large to small. We have not 
discussed in this chapter the links to solar neutrinos nor the LSND 
results pot. Certainly there is no conflict between the atmospheric muon neu- 
trino results and the possible (nay likely) solar oscillations. If, however, both 
solar and the LSND results are correct, then we have surely some inter- 
esting physics to untangle, as it is generally admitted that no simple three 
neutrino model can incorporate all three neutrino anomalies, and that new 
degrees of freedom would be required (see Chapter IX). From the evidence 
presented in this chapter alone, however, we can make some far reaching, 
perhaps paradigm shifting, conclusions. 

1.3.1 Astrophysics and Cosmology 

The implications of the oscillations results have been explored in other chap- 
ters in this book, so here we only outline those relative to the muon neutrino 
oscillations. First, it appears that neutrinos with summed masses of the or- 
der of 0.1 eV will not make any major contribution to resolving the dark 
matter quandary. Nonetheless with a ratio of 2 billion to one for photons 
(and neutrinos) to nucleons from the Big Bang, even such a small neutrino 
mass may be greater in total than all the visible stars in the sky. Hence while 
one must account for neutrino mass in further cosmological modeling, neu- 
trinos are not likely to constitute the bulk of the "missing matter" . However 
if the neutrinos should be nearly degenerate in mass and all have masses 
in the range near 1 eV (and hence we are observing only small splittings 
with the oscillations), then neutrino mass may dominate the universe. While 
neutrinos are not favored by astrophysical modelers (fitting the spatial fluc- 
tuations in the cosmic microwave background for example), large neutrino 
masses are not ruled out (SMi < 6 eV, Ho = 65 km/s/Mpc, i = 1, 2, 3)||5]|. 



1 Atmospheric Neutrino Anomaly 33 

Nearly degenerate neutrino masses would not present a consistent picture 
with the quark and charged lepton masses, which make large mass jumps 
between generations. But who knows? We do not have a viable GUT with 
mass predictions so an open mind is appropriate. 

The other major area of significance, perhaps of the deepest significance, 
has to do with baryogenesis, the origin of the predominance of matter by one 
part per billion over anti-matter at Big Bang time. There are claims that the 
old idea of accumulation of net baryon number via CP violation in the quark 
sector, while satisfying the Sakharov conditions |^3| , may not sufiice from the 
early stages of universe expansion |]5^. If that is indeed the case, it may be 
that neutrinos provide the avenue for net baryon asymmetry generation, with 
the expression into hadrons becoming manifest relatively late in the game at 
electroweak phase transition ||55||. CP and even CPT violation in the neutrino 
sector, as yet almost unconstrained, could have dramatic implications. 

Neutrino masses and possible sterile neutrinos have also been invoked to 
help resolve problems in understanding heavy element synthesis in supernovae| 



1.3.2 Theoretical Situation: Why So Important 

Other chapters in this volume deal with the particle theory situation, so we 
make only several general remarks here, more from the experimentalists' phe- 



nomenological viewpoint. Figure 1.24, shows the masses of the fundamental 



fermions in three generations on a logarithmic scale in mass. Dramatically, 
one sees that if the neutrino masses are near the lower bounds (that is at the 
presumed mass differences from present atmospheric and solar results) , they 
lie 10-15 orders of magnitude below the other fundamental fermions (charged 
quarks and leptons). Graphically one notes the spacing between the neutrino 
masses and the charged fcrmion masses is just about the same as the distance 
(on the log scale) to the unification scale. This is a pictorial representation 
of the see-saw mechanism, as noted more than ten years ago|Q. This points 
up the task for grand unification model building, and highlights the deep link 
between neutrino masses and nucleon decay [p2|. 



1.3.3 Future Muon Neutrino Experiments 

The physics community seems to have rapidly accepted, with caution, the 
seeming inevitability of neutrino mass and oscillations [Q. Yet most proba- 
bly the game has hardly begun and many a subtlety may await our explo- 
ration. However if the LSND claims will go quietly away (after the BOONE 
experiment runs), the mass and mixing picture could settle into the simple 
hierarchical pattern as explored in the bi-maximal mixing scenario (or similar 
versions). This highlights the importance of experiments to follow up on the 
LSND results as one of the first agenda items in the current neutrino business 
(see Chapter VII). 
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Fig. 1.24. The masses of the fundamental fermions. The lower shaded region in- 
dicates the rough range allowed by present oscillations results, with the upper 
boundary for nearly degenerate neutrino masses and lower bound for approximate 
atmospheric and solar results with an extrapolation to the first generation mass. 
The unification scale marks the rough range implied by the see-saw mechanism. 



Given present indications, it would seem that the K2K|5^, MINOS[|60{ 
and CERN-Gran-Sasso experiments ICANOE[|l) and 0PERA(62| should 
confirm the SuperK results and make the oscillation parameters more precise. 
Of course, one would really like to see tau appearance, not just muon disap- 
pearance to be sure we are not being misled. There has been some debate in 
the community as to what constitutes appearance. Because of the complex- 
ity of tau final state identification, this author would prefer to see a real tau 
track recorded. In any case plans are in progress in the US, Japan and Europe 
for the obvious follow up experiments to solidify the oscillation scenario and 
refine the parameters, which experiments should take less than a decade. At 
this time it appears that only OPERA has the opportunity to record physi- 
cal tau lepton tracks in emulsions, though the other experiments (including 
SuperK) may be able to detect tau kinematic signatures statistically. 

More interesting for the long range physics is filling in the MNS ma- 
trix (lepton equivalent of quark CKM matrix) for neutrinos This is not 
an easy business. The atmospheric neutrino measurements really are only 
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defining, at best, three of the nine elements! Solar neutrinos get us another, 
perhaps a constraint on two. Measuring the tau related components directly 
seems pretty hopeless. Of course if we can assume the matrix to be unitary 
and real we would be in good shape as there are then only three independent 
parameters (plus the masses). But we do not know this, and if there exist CP 
violations we then have a total of three angles and one phase plus possibly 
another phase for the case of Majorana neutrinos, but which phase is only 
measurable in special circumstances such as neutrino-less double-beta-decay. 
If there are more (heavy or sterile) neutrinos, then things could be much 
more complicated (as the 3 by 3 sub-matrix will not be unitary) . By analogy 
with the quarks (where the 3 by 3 CKM matrix with small mixing angles 
and one CP violating phase suffices), perhaps we should not worry too much, 
except for lack of any guidance whatsoever from theory. CP violation is only 
very weakly constrained experimentally in the neutrino sector at present, 
so we could be in for big surprises, and given the neutrino connection with 
cosmology and baryogenesis, one should indeed be suspicious, I believe. 

As a whole, the particle physics community is just beginning to examine 
the newly illuminated possibilities in the neutrino sector. At the moment 
it appears as though muon colliders may provide our best route for next 
generation explorations in this realm p^ , |65[ . The serendipitous realization 
that muon storage rings can provide neutrino beams more intense by about 
4 orders of magnitude than previous artificial neutrino sources, along with 
the potential, via polarization of these beams, to make nearly pure beams 
of muon neutrinos or anti-neutrinos, with a relatively narrow energy spread, 
allows one to dream of experiments not heretofore thought possible. The 
prospects for this endeavor are just emerging at the time of this writing and 
the dialog between physics possibilities and machine technical capabilities is 
in much ferment. The present excitement about this endeavor, probably at 
least ten years from realization, promises much evolution of thinking about 
critical tests and experiment planning. Aside from simple checking of previous 
results, precision measurement of oscillation parameters, and filling out of the 
MNS matrix, the major goal in this authors mind is exploration for CP and 
CPT violation in the neutrino sector, not possible by any other means yet 
conceived. Of course if ^13 turns out to be zero, which hopefully we will know 
in a few years, then CP violation is not be possible. Credibly motivated CPT 
violation can still be a possibility however |]67|. 



Table 1.2 shows a rough comparison of the physics capabilities and reach 
of SuperK and the various present and future long-base-line experiments. 
This table was generated by an informal working group at the 2000 Aspen 
"Neutrinos with Mass" Workshop. Although the details may annoy propo- 
nents of some projects (and please others), the purpose was not to evaluate 
specific experiments, but to try and gauge the strengths of the various ap- 
proaches as we move forward in the quest to complete the MNS matrix and 
explore for CP and CPT violations in the neutrino sector. We had inadequate 
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information about OPERA and the proposed 50 GeV proton driven neutrino 
beam in Japan (JHF to SuperK), so they may not be fairly represented. 

The table assumes that the Am"^ = 0.003 eV'^ and sin^(26') ~ 1.0 with 
dominantly v^, ^ i>r mixing, large mixing angle MSW solar neutrino oscil- 
lations and no LSND/BooNE indication of oscillations. We also assumed a 
20 — 50 GeV neutrino factory and high quality neutrino detector at a few 
thousand kilometers distance. 

There was considerable discussion of whether the new experiments will 
or will not be able to see at least one dip in the detected neutrino spectrum, 
unambiguously discriminating oscillations from disappearance (as in decay). 
My conclusion is that such a detection will be very difficult for all experiments 
with traditional neutrino beams if the Arri^ is at the low end of the allowed 
region. ICANOE claims to be able to detect oscillations in the cosmic ray 
beam, a consequence of their excellent resolution and the larger range oi L/E 
values available from the atmospheric neutrinos than from a fixed distance 
long-base-line beam. 

I think what we all learned from this exercise is that there really is a great 
gulf between what we can accomplish with a neutrino factory and anything 
prior to that, even with more powerful traditional neutrino beams. 



Table 1.2. Comparison of SuperK and long-base- line neutrino experiments in terms 
of physics capabilities and reach. See text for explanation. 



Experiment: 


SuperK 
(atm v) 


K2K 


MINOS 


ICANOE 


OPERA 


JHF2K 


V Fac 


d{ln{Am^)) 
d(sin^ 261) 


50% 


20% 


10% 


10% 


20% 


10% 


2% 


5% 


5% 


5% 


5% 


? 


4% 


2% 


See Oscillations? 


X 


X 


v/(?) 


X (atm^O 


? 


X 




r appear (kink) 


X 


X 


X 


X 




X 




T appear (kinem) 


^? 


X 


V 


3V 




X 




sin^ 2013 Umit 


0.1 


0.03? 


0.03 


0.015 


(ident e?) 


0.03 1 


- 3 X 10 


d{vs/nr) 


20% 


X 


5%? 


5% 


? 


X 


1% 


Elim Decay Models? 


x(? NC/CC) 


X 


v/(NC/CC) 


X (atmi/) 


? 


X 




Sign of Arn^ 


X 


X 


X 


X 


X 


X 




VfL Vt 


X 


X 


X 


X 


X 


X 


^/ 


CP violation tests 


X 


X 


X 


X 


X 


X 




CPT violation tests 


X 


x(?) 


7 




X 


x(?) 





Measuring absolute neutrino mass remains a frustrating problem, which 
will not be resolved in the near future it seems. While pushing to lower mass 
limits with Tritium beta decay experiments will apparently not be able to 
reach below 0.1 eV ^ there is some hope from CMBR measurements JsTf , there 
is a long shot via the "Weiler process" , and optimistically neutrino-less 
double-beta-decay experiments may eventually reach 0.01 or even 0.001 eV . 
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At high energies, explorations for cosmic neutrinos may be carried out 
in deep arrays in the ocean and under icej6|]. While the main goal of these 
attempts at high energy neutrino astronomy will be aimed at astrophysics, 
with a high-energy-threshold-detector capable of registering neutrinos in the 
PeV range, it may be that such instruments will be able to directly detect 
tau neutrinos (via the "double bang" signature |]69{), and even determine the 
neutrino flavor mix, to the benefit of both particle physics and astrophysics. 

A next generation (megaton) scale nucleon decay instrument to probe 
lifetimes to 10"^^ years would do wonders for advancing neutrino physics as 
well. Simply building a larger version of SuperK will not suffice because of 
the need for greater resolution as well as size. The only candidate I see to 
go beyond SuperK is something like the AQUA- Rich style of imaging water 
Cherenkov detector js^,^. An attractive alternative which need not be so 
massive to get to 10^^ years in kaon modes of nucleon decay, might be a 
50 — 70 kiloton liquid argon detector of the ICARUS style. Perhaps such a 
detector can be realized in concert with a long baseline beam from a neutrino 
factory. 

From the foregoing it should be apparent that we have entered a new era 
in elementary particle physics, and that one can expect a long and interest- 
ing exploration into neutrino mass and mixing now that the door has been 
opened. 
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